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ABSTRACT 

Recent studies based on line indices of nearby early-type galaxies suggest that stars in the galaxies' 
inner part may be younger and more metal rich than stars in the outer part. If confirmed, the finding has 
a profound implication for the formation and evolution of early-type galaxies, since such age gradients 
naturally arise if merging played an important role in their formation/evolution. As an indpendent 
test for the existence of the age gradient, we investigate the optical-NIR color gradients of six field 
early-type galaxies at z ~0.4 - 1.0 where the age gradient is expected to be detectable when combined 
with z = calibrations. By utilizing HST NICMOS H-band imaging along with WFPC2 optical data, 
we obtain a broad wavelength baseline, giving results that are sensitive to the expected age gradient. 
When compared with simple model predictions, five of the objects show negligible age gradients and 
are metallicity gradient dominated, while the remaining object shows a color gradient consistent with 
the age-|-metallicity gradient picture. Our results suggest that stars within early-type galaxies may be 
coeval at z ~ 1, but sprinkling of young stars might have occured at z < 1. 

Subject headings: cosmology: observations - galaxies: elliptical and lenticular, cD - galaxies: evolution 



look as blue as, or bluer than, the outer region if the 
age-|-metallicity gradient model is correct. Consequently, 
the color gradient would become flatter or even reversed 
as a function of rcdshift. 

Color gradients of distant early-types have been studied 
using optical HST images, and such studies have excluded 
models where the color gradient is driven purely by age 
(Tamura et al. 2001; Saglia et al. 2000). However, they 
put few constraints on the more plausible age-|-metallicity 
gradient model, due to the fact that the age effect, when 
buried under the metallicity gradient, is barely detectable 
in observed optical passbands (see section 5). Several stud- 
ies find E/SO candidates with bluer cores (e.g., Abraham 
et al. 1999; Menanteau et al. 2001; Im et al. 2001a), but 
a detailed study of their kinematic and structural proper- 
ties shows that these galaxies do not appear to evolve into 
typical E/SOs today (Im et al. 2001b). 

In an attempt to test the age-|-metallicity gradient 
model, we have studied optical-NIR (near-infrared) color 
gradients of early- type galaxies out to z 1, using HST 
NICMOS data for NIR and HST WFPC2 data for optical. 
In this paper, we present our preliminary results which 
favor no significant age-gradient in the majority of early- 
type galaxies in our sample. 

2. SAMPLE 

Our sample consists of six field E/SO galaxies at z < 1 
of which five are located in the Groth strip (Rhodes et 
al. 2001), and one from the Hawaii Deep Field (Cowie et 
al. 1994). These objects come from the NICMOS snap- 
shot survey of distant galaxies (HST G07895 program). 
The snapshot survey obtained HST NICMOS Camera 2 
iJ-band data (F160W, exposure time of 1200 sees per ob- 
ject) of 16 field galaxies at 0.2 < z < 1 for which redshifts 

^ Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
^ Physics Department, University of California, Santa Cruz, CA 95064; shinkley@ucolick.org 
Astronomy & Astrophysics Department, University of California, Santa Cruz, CA 95064; myung@ucolick.org 
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1. INTRODUCTION 

Early-type galaxies in the local universe show a mod- 
erate change in stellar populations from their centers to- 
ward their edges. Stars are redder toward the center of 
early-types (e.g., Franx & lUingworth 1990; Peletier et al. 
1990a, 1990b; Kormendy & Djorgovski 1989), and stud- 
ies of metallicity-sensitive line indices of nearby ellipticals 
(e.g., Kobayashi & Arimoto 1999; Gonzalez 1993; Car- 
ollo et al. 1993) show that a change in metallicity has a 
substantial effect on the color gradient. However, when 
age-sensitive Balmer line indices (see Worthey 1994) are 
investigated, some ellipticals not only show metallicity gra- 
dients, but also age gradients in which the inner part is 
younger than the outer part of the galaxy (Gonzalez 1993; 
Tantalo, Chiosi, & Bressan 1998; Trager et al. 2000). 

This picture of an age-|-metallicity gradient makes sense 
if mergers play a significant role in the formation/evolution 
of early-type galaxies. Only the metallicity gradient 
is expected in models where early-type galaxies formed 
via monolithic collapse at high redshift and evolved pas- 
sively thereafter (e.g., Carlberg 1984). The indication of 
younger stars near the center suggests additional star for- 
mation/accretion activities very possibly through merging. 
However, a number of studies point out difficulties of using 
Balmer lines for age measurements (Maraston & Thomas 
2000; Bressan et al. 1996; Lee, Yoon, & Lee 2001; Jor- 
gensen 1997; Fisher, Franx, & lUingworth 1995) — thus, 
an independent assessment of the age-metallicity gradient 
model is highly desired. 

One promising way to detect age gradients is to study 
color gradients of distant early-type galaxies (z > 0.2). 
As one looks back in time, younger stars become more 
visible out of old stars where they are embedded. There- 
fore, the inner parts of distant early-type galaxies would 
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as well as HST optical data are available in F(F606W) 
and/or /(F814W) bands. 

E/SO galaxies are chosen using the quantitative selection 
criteria of Im et al. (2001a; 2001b) imposed on the /-band 
data. The method uses two parameters: bulge-to-total 
light ratio (B/T) which measures the significance of the 
bulge component, and the residual parameter (R) which 
measures the smoothness and symmetry of the galaxy. We 
adopt B/T > 0.4 and R < 0.08 following Im et al. (2001a). 
Note that these parameters are derived for each galaxy in 
the Groth strip using the 2-dimensional surface brightness 
fitting algoritm, GIM2D (Simard et al. 2001), which also 
returns other parameters such as the half light radius (rhi). 
For one galaxy in the Hawaii Deep Field, the B/T value is 
taken from the MDS database (Ratnatunga, Ostrander, & 
Griffiths 1999). Although an R value is not readily avail- 
able for this object, a visual inspection confirms that it is 
an E/SO. The basic data for our sample are presented in 
Table 1, and images are shown in Fig. 1. Note that magni- 
tudes are in the Vega system for the HST filters (Holtzman 
et al. 1995; Simard et al. 2001). 

3. DATA REDUCTION AND ANALYSIS 

The reduction of the iJ-band data was performed using 
the IRAF'' NICPROTO package (Bushouse et al. 2000). 
The NICPROTO package^ includes the standard NIC- 
MOS pipeline procedure (calnica and calnicb), but it gives 
special attention to removing image anomolies such as 
a pedestal variation across the array and a varying bias 
senstivity in the detector itself. In addition, all of the 
data were reduced using the new NICMOS synthetic dark 
frames, which helped to alleviate a time-dependent shad- 
ing feature caused by the readout amplifiers (Monroe & 
Bergeron 1999). We have carried out these additional 
reduction procedures in order to remove the background 
variation as much as possible, since color gradient mea- 
surements require very carehil background determination. 

Four dither exposures were obtained for each field. From 
these four, a single higher signal-to-noise image was pro- 
duced after the four were cross-correlated and shifted using 
subpixel translations with cubic spline interpolation. Sim- 
ulations using tinytim PSFs confirm that the simple shift- 
and-add method used here will not substantially change 
the characteristics of the NICMOS PSF. The'WFPC2 
images were obtained from either the DEEP or MDS 
databases. Once the calibration and reductions were com- 
plete, the H-band Vega magnitudes were calculated using 
the relation H = 21.75 - 2.5 log(i:' TV/sec) (Schultz et al. 
2001), where DN (counts) is obtained from GIM2D. 

Since NICMOS and WFPC2 have different PSFs, we 
matched the PSF by cross-convolving our data - i.e., the 
WFPC2 images were convolved with the NICMOS PSF 
and vice versa. The PSFs were created at the object lo- 
cation using the TinyTim package (Krist 1993). Since the 
plate scales of the WFPC2 (0'.'0996)and NICMOS Cam- 
era2 (0'.'075) are different, oversampled PSFs were pro- 
duced and then rebinned to have the platescale of the 
other detector. Note however, that no cross-convolution 
was performed for the {V — I) gradients since the differ- 
ence in the V and I PSF is insignificant. 

^ IRAF is distributed by the National Optical Astronomy Observatories, under contract to AURA, Inc. 
5 The most up-to-date NICMOS tasks in STSDAS now include most NICPROTO procedures. 



The standard IRAF photometry package "ellipse" was 
used to fit elliptical isophotes to each object, and care 
was taken to ensure that the ellipticity and position an- 
gles were consistent between the NICMOS and WFPC2 
data. Moreover, each object's isophotal centroid agreed 
extremely well between the two bands, and these were con- 
sistent with the GIM2D values to within 0.3 pixels. 

The background value was determined by taking a me- 
dian sky value of the region beyond 5 rhi and within a 6'.'4 
by 6'.'4 square region of the object. We checked the sky 
value by constructing the flux growth curve, and found 
that the growth curve nearly flattened out beyond 4 ru. 
The sky value is consistent with the output from GIM2D 
which uses a similar procedure (<10%). Slight differences 
in sky values affect the aperture photometry at the outer 
part of the galaxy at 3 - 5 r^j;, thus we have restricted our 
study to regions inside r < 2rhi. 

4. MODELS 

This study uses simplifled versions of the two kinds of 
models discussed in section 1. The first is a pure metallic- 
ity gradient model, and the second is an age-|-metallicity 
gradient model. The metallicity {Z) and the age (t) gra- 
dients are modelled as: 

Slog{Z) = C + Czlog{r/rhi), (1) 
6\og{t) = C' + Ct\og{r/rhi), (2) 
where C and C are constants, rhi is the half light radius, 
and Cz and Ct are the slope of the metallicity gradient 
and the age gradient respectively. 

The pure metallicity gradient model assumes no age gra- 
dient (i.e., Ct = 0), and local E/SOs show Cz values in the 
range of -0.4 < Cz < -0.1 with (Cz) = -0.2 (Henry & 
Worthey 1999; Kobayashi & Arimoto 1999). The plausible 
range for the central metallicity value is 100%-300% 
within rhi/8 (Trager et al. 2001), and we use 250% for 
the inner region defined at log(^) ~ —1.175 (the light- 
weighted mean radius within rhi/S). The outer region is 
defined as where the metallicity falls to 100% Z^. The 
model colors are calculated at these inner and outer radii, 
and the model color gradient is the straight line which con- 
nects the two points. Then, we find the best-fit Cz and 
the age (or Zfor, the formation redshift) by minimizing a 
fit to the observed color gradient. 

The age-|-metallicity gradient model assumes an age gra- 
dient in addition to the metallicity gradient. Line indices 
studies find that (Ct) ~ 0.1 and '{Cz) ^ -0.25 for local 
E/SOs (Henry & Worthey 1999). Our base age-|-metallicity 
gradient model assumes these slopes, and the zero point 
of the base model is adjusted to fit the optical-NIR color 
at r/rhi ^0.1 0.2 by changing Zfor- Besides the base 
model, we allow our model to fit for Ct, Cz and Zfor- 
Again, 250% Z^ is placed in consideration for the inner 
metallicity. Note that we use only two points for the color 
gradient predictions. Consequently, the predicted color 
gradients are straight lines. This should be a good first or- 
der approximation for the color gradients given the quality 
of our data, although the color gradients can get curved 
when the age effect is strong (e.g., Tamura et al. 2001). 
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GSS: ()52_6543 z=0.76 



SSA13 55 z=1.03 GSS: ()82_6533 z=(>.462 




Fig. 1. — The /-band (top row) and the //-band (bottom row) images of E/SOs in our sample. All images are oriented the same way and 
are on the same scale. 



SEDs arc generated at different ages and nietallicities 
using the 1996 version of Bruzual & Chariot model (1993). 
The standard Salpeter IMF, 0.1 Gyr burst, and an open 
universe with Qq = 0.2 and Hq = 70 km sec~^ are as- 
sumed. The SEDs are redshifted and Jr— corrections are 
derived using the formula in Im et al. (1995), which are 
then used for calculating colors. Note that the zeropoint 
of the optical-NIR color may vary with spectral synthe- 
sis models, but the slope of the color gradient is not likely 
to be model-sensitive (Chariot, Worthey, & Bressan 1996). 
Also, the predicted slope of the color gradient is insensitive 
to the assumed cosmology. 

5. RESULTS 

Color gradients for our sample are shown in Fig. 2. Ex- 
cept for one object without F-band data, two plots are 
shown for each object — one for the V — H gradient and 
another for the V — I gradient. At 2; = 1, the observed 
V—H gradient corresponds roughly to the rest-frame U—I 
gradient. Also plotted in Fig. 2 are model predictions. 

The first five objects have V — H 01 1 — H gradients in- 
consistent with that of the base age-|-metallicity gradient 
model predictions, while their color gradients can be well 
matched by the pure metallicity gradient model with the 

best-fit slope of Cz = (-0.2 0.4) ± 0.05. Fig. 2 also 

demonstrates that the V — I gradient is not very sensitive 
to the age gradient, thus not useful for testing the age- 
metallicity gradient model. Since these objects are already 
at the lookback time of ~ 7 Gyrs, the mean stellar age in 
the inner part must be less than a few Gyrs under the 
base age-hmetallicity gradient model. This would produce 
flat or reversed optical-NIR color gradients (dashed lines) , 
which does not match the data. Thus, it is unlikely that 
any significant secondary star formation/accretion event 
occured at the center of these galaxies before z ~ 1. To 
prove this point further, we fit the color gradient for Ct 
using a fixed Cz = —0.25, and find no significant age 
gradient (Ct ~ 0.01 ± 0.02). One might consider a case 
where the age gradient is hidden under a steep metallic- 
ity gradient. To check such possibility, we fit for Ct using 
Cz = -0.5 and find that C( ~ 0.08 ± 0.02 for the five 
objects. The derived Ct values are inconsistent with the 
local constraints where Ct — 0.3 when Cz — —0.5 (Henry 
& Worthey 1999). Thus, we still require more young stars 
near the center for this picture to be true. Note that the 



model color is sensitive to the assumed metallicity. If one 
uses a different metallicity value, the zero point of the 
predicted color gradient shifts too. However, this shift in 
the zero-point can be compensated by changing Zfor- We 
tried 100% for the metallicity of the central part. This 
changed the predicted color gradient slightly, but did not 
affect the overall conclusion about the age gradient. 

The {V — H) color gradient of the last object 
(GSS082_6533) shows that its inner part is bluer than 
the outer part, with the best-fit gradients being Cz = 
-0.27 ± 0.20 and Cj = 0.20 ± 0.10. This is consistent 
with the age-hmetallicity gradient picture, and sprinkling 
of young, blue stars might have occurred in the past few 
Gyrs, near the center of this galaxy. 

Fig. 2 also shows that the sensitivity to the age gra- 
dient seems to be; higher for E/SOs at high-z. For ex- 
ample, the difference between the pure-Z model and the 
age-|-metallicity gradient models is more drastic for the 
z = 0.89 object than the z = 0.462 object. However, sec- 
ondary star formation may occur after 2: ~ 1, in which 
case it is advantageous to observe galaxies at moderate z. 

6. CONCLUSIONS 

We have performed surface photometry on six early-type 
galaxies at z = 0.4 1.0 and obtained their optical-NIR 
color gradients. These were then compared with simple 
models driven purely by metallicity or models driven by 
changes in age and metallicity. The age-|-metallicity model 
predicts a flattened or reversed optical-NIR color gradient 
due to the assumed existence of younger stars in the inner 
region. However, the model fails to reproduce the observed 
color gradient of the majority of our galaxies (5/6). There- 
fore, we conclude that no significant secondary star forma- 
tion/accretion occured at 2: > 1 for these galaxies. While 
evidence is lacking for a strong age-gradient for z ^ 1 
early-types, one object at the lowest redshift {z = 0.462) 
shows a color gradient consistent with the age-|-metallicity 
gradient model. This shows that the age-gradient is de- 
tectable, and that secondary bursts/accretions may occur 
at moderate redshifts. In the future, we hope to draw 
firmer conclusions on this issue by enlarging our sample 
and improving our model predictions. 

This work was supported by the STScI grant GO- 
07895. 02-96A, and partial support was provided by the 
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GSS 044_1676, z=0.90 
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Fig. 2. — Shown here arc pairs of {V — H) and (V — I) gradients for all of our data. Object SSA13 #55 docs not have observations in V, so 
only ail I — H gradient is shown. The solid lines indicate the purely metallicity-driven model, while the dashed lines are the age -|- metallicity 
models. Also shown in the figure are the best- fit or assumed Cz, C't, and the age for the models. For 052_6543, a slightly shallower age 
gradient is used (0.05) since the canonical value produces a color gradient that is negatively too steep. For a majority of this sample (5/6), 
the gradients follow the purely metallicty-driven model. Only the last object, GSS 082_6533, shows a sign of young stars in its inner region. 
In all of the plots, the points are closely spaced, so two adjacent isophotes may share some pixels. Thus, the error bars arc not independent 
of each other. The lookback times of 6, 7, 8, 9, and 10 Gyrs correspond to redshifts of 0.79, 1.1, 1.5, 2.2, and 3.4 respectively. 
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Table 1 

Basic Data for Early- type sample 



Object ID 


RA 


DEC 




Type 


H 


I 


V 


(B/T), 


(10) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


O) 


(11) 



GSSa52_6543 

GSS044_1570 

GSS044_1676 

GSS053_6418 

SSA13 #55 

GSS082_6533 



14:17:48.2 
14:18:04.8 
14:18:04.9 
14:17:55.6 
13:12:19.3 
14:17:28.8 



52:31:17.3 
52:30:59.1 
52:30:53.6 
52:31:52.4 
42:45:01.2 
52:27:37.8 



0.756 
0.89 ± 0.15* 
0.90 ± 0.15* 

0.675 
1.03 

0.462 



SO 
E/SO 
E/SO 
Sp.SO 

SO 
E/SO 



18.02+;0| 

18.96j:i| 
19.82+-10 

18.53l;Jl 
18.45±;1« 



20.31±-21 

2i.76±:ga 

+ .05 



22.29 
20.97 
21.82 
20.57 



05 
+ .01 
-.02 
+ .22 
.22 
+ .01 
01 



22.26+'^| 
23.79±:°^ 
24.33+-5q 
22.66^°^ 

N/A 
22.04+;°! 



44+-°l 

71+ oe 

' — .06 

o.56+;i| 

72+ °=* 
0.63111 
41+ °^ 



0.637+-^2 
0.175i:°^^0 
0.195+0=1 
0-354i:°l| 

0.379+;°!^, 

0.207±:°J° 



0.07 
0.06 
0.04 
0.04 
N/A 
0.05 



Note. — (1) Source ID given by FFC-XXYY, where FF is the Groth Strip subficld, C is the WFPC2 chip number, and XX and YY arc the chip 
coordinates in units of 10 pixels. The "SSA" denotes that the object is from the Hawaii Deep Field. (2) Right Ascension (J2000). (3) Declination 
(J2000) (4) Spectroscopic redshifts(5) IMorphological type (6)-(8) Total magnitudes for H, I, and V-bands. (9)- ( 1 l)BuIge-to-totaI ratio, half-light 
radius, and residual parameter within 2 rf^i for the /-band data as me;isured by GIM2D (Simard et al. 2001 in preparation). 



* Photometric redshifts are derived as described in Im et al. (2001a) 



